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I. INTRODUCTION 



The process studied in this paper consists of two steps. First, a free electron is recombined 
into an electronic shell of a positive, preferably highly charged ion while resonantly transferring 
its excess energy to the nucleus. This step is referred to as nuclear excitation by electron capture 
(NEEC) in the literature. It is the time-reversed process of internal conversion (IC) and can also 
be regarded as the nuclear analogue of dielectronic recombination (DR), where a bound electron 
is resonantly excited. In the second step, the excited nuclear state thus formed decays radiatively. 
The whole process is illustrated in Fig. \T\ Proposed for the first time in Ref. [1], the NEEC 
recombination mechanism has not been observed experimentally yet. 

Related processes at the borderline between nuclear and atomic physics have been experimen- 
tally confirmed, such as nuclear excitation by electronic transition (NEET), that is a decay mode 
of the electron atomic shell in which energy is transferred to the nucleus. NEET has been observed 
in 189 0s [2], 237 Np |3], and, most recently, in 197 Au [4]. For its time-reversed process, the bound 
internal conversion (BIC), a direct experimental evidence was not found until recently qO. 

The detection of the recombined ions or the emitted radiation characteristic for the nuclear 
decay can be used to observe NEEC experimentally. The position of the peak in the recombination 
cross section or the spectrum of emitted photons directly yields the nuclear transition energy. 
The detected spectra may provide further information on nuclear decay mechanisms and level 
population. Thus, the study of NEEC in recombination experiments can also be regarded as a 
possible novel tool to investigate nuclear properties. 

In a previous paper flg], we presented total NEEC cross sections for a range of heavy isotopes. 
Recently, the quantum interference of NEEC and the background process of radiative recombina- 
tion (RR) was investigated [7fl. In the present work we study a new aspect, namely, the alignment, 
i.e., the magnetic sublevel population of the nuclear excited states generated in the resonant re- 
combination process. The alignment of the nuclear excited state also allows one to derive the 
angular distribution of the radiation emitted in the nuclear relaxation process and the cross section 
differential with respect to the photon emission angle. Differential cross sections for the similar 



process of DR have been investigated theoretically and experimentally 11121. [ 

The photon angular distribution is of general interest for experimental implementations aiming 
at observing NEEC by the detection of radiation. Furthermore, as found in this work, the photon 
emission pattern is substantially different from the one of the concurrent process of RR, which 
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FIG. 1: NEEC recombination mechanism of a continuum electron into the K shell of a bare ion. The 
nucleus is schematically represented as undergoing the transition from the ground state (G) to the excited 
state (E) and again to its ground state. The figure below shows the geometry of the process. 

may help to discern the two channels. 

The reorientation of the nuclear axis caused by the electric field of a charged particle that excites 
the nucleus is a well-studied process in nuclear physics. The change in the nuclear spin directions 

lisions with low-energy charged particles 



following the Coulomb excitation of the nucleus by co 
affects the angular distribution of the emitted 7 rays II 1411 . In the case of the electric radiative 
transitions, NEEC can be regarded as a Coulomb excitation with free electrons that are in the end 
recombined into a bound state of the ion. 

We apply a density matrix formalism to describe the alignment of rotationally excited nuclear 
states formed by NEEC and the angular distribution of the de-excitation radiation. Nuclear states 
are treated in the framework of the nuclear collective model and the electrons are described by 
Dirac four-component wave functions, as necessary for the heavy elements studied here. The 
electron-nucleus interaction operator is expanded in spherical multipoles to facilitate the evalu- 
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ation of transition matrix elements. This theoretical approach is presented in detail in Sec. HQ 
NEEC nuclear alignment parameters are presented for the case of E2 excitations in heavy highly 
charged ions and differential cross sections are compared to the same quantity corresponding to 
the background process of RR in Sec. [TTTJ The relevance of these results for future experiments 
targeting the observation of NEEC is discussed. We briefly summarize the findings of the paper in 
Sec.|IVl Atomic units are used throughout this work unless otherwise specified. 



II. THEORETICAL FORMALISM 



Since its introduction in 1927 by von Neumann and Landau, the density matrix approach has 
been found to be a useful and elegant tool in many fields of modern physics. For applications in 
atomic physics, and combined particularly with the concept of spherical tensors, this approach was 
originally developed by Fano jl5\\ in the 1950s. Since then, the density matrix theory has been 
utilized successfully in many case studies on atomic collisions, e.g., for describing the excitation of 
atomic autoionizing states, polarization effects in radiative and Auger decays, cascade processes, 
or lifetime interferences in resonantly excited atoms. 

In this paper we apply a density matrix formalism to describe the two-step process in which 
(i) a free (or quasi-free) electron is captured into the bound state of an initially bare ion with 
excitation of the atomic nucleus which (ii) subsequently decays under the emission of character- 
istic radiation. Because the properties of this radiation are closely related to the alignment of the 
excited nuclear state, we first have to investigate the population of these states as it arises due to 
the electron capture process. Therefore, in Section ITTAl we derive the general formulas for the 
density matrix of the excited nuclear states produced by electron capture. In particular, we show 
how the sublevel population of these states can be described in terms of the so-called alignment 
parameters. The calculation of these parameters involves the NEEC transition amplitudes and re- 
quires the use of a nuclear model. Following Ref. [6], the NEEC transition amplitudes are derived 
in Section MB I Finally, in Section III CI we consider the subsequent nuclear decay and obtain the 
angular distribution of the de-excitation photons with the help of the alignment parameters. 
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A. Alignment of the excited nuclear state 

Within the density matrix theory, the state of a physical system is described in terms of sta- 
tistical (or density) operators. These operators can be considered to represent, for instance, an 
ensemble of systems which are — altogether — in either a pure quantum state or in a mixture of dif- 
ferent states with any degree of coherence. Then, the basic idea of the density matrix formalism is 
to accompany such an ensemble through the collision process, starting from a well-defined initial 
state and by passing through one or, possibly, several intermediate states until the final state of the 
collision process is attained. 

In NEEC, the initial state of the combined system is given by the electron with a well-defined 
asymptotic momentum p and spin projection m s , and an ion which is specified in term of its 
nuclear charge Z and its initial nuclear spin quantum numbers 7, and projection Mj. (We note 
here that we consider recombination into bare or closed shell ions, thus the angular momentum 
of the bound electron shell is zero.) Assuming that these two subsystems, ion and electron, are 
uncorrelated, the overall initial density operator is given as the direct product of the two initial 
subsystems' density operators: 



If neither the electrons nor the ions are polarized initially, an averaging over the magnetic quantum 
numbers can be performed and the tensor product can be written as 



where iV; denotes all the additional quantum numbers needed for a unique specification of the 
nuclear states. 

In the intermediate state d formed by the capture of the electron, the statistical operators have 
to describe both the electron in some bound ionic state \n d n d j d m d ) as well as the state of the 
excited nucleus \N^I d M d ). Here, n d , n d , j d and m d are the principal quantum number, Dirac 
angular momentum quantum number, total angular momentum quantum number and magnetic 
quantum number of the bound one-electron state, respectively. As known from density matrix 
theory, the statistical operators of the initial and the (subsequent) intermediate states of the system 
are connected by 



Pi = Pion ® Pe ■ 



(1) 




Pd = T en PiTl 



0) 
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where T en is the transition operator for the electron-nucleus interaction which causes the nuclear 
excitation. The particular form of this operator will be given in Section ITTBl 

Instead of applying equation ©, in practice, it is often more convenient to rewrite the statistical 
operators in a matrix representation. For instance, in a basis with well-defined angular momenta, 
the intermediate- state density matrix is given by 

(N d I d M d , n d n d j d m d \p d \ N d I d M' d , n d n d j d m f d ) = 

~ 2j~^i ^2 { N *d!dM d , n d K d j d m d \T en \ NJiMi, pm s ) 

1 m a Mi 

x (N*I d M' d , n d K d j d m' d \T en \ NJiMi, pm s )* , (4) 



assuming that both the incident electrons and ions are initially unpolarized [see Eq. ©]. Indeed, 
the intermediate-state density matrix © still contains the complete information about the NEEC 
process and, thus, can be used to derive all the properties of the bound electron and the excited 
nucleus. For instance, assuming that the magnetic states m d of the bound electron remain unob- 
served in the particular experiment, we may characterize the sublevel population of the excited 
nucleus \N d I d ) in terms of the nuclear density matrix 

(N* d I d M d \ P T d \N* d I d M' d ) = 2 (NZIdM^niKdjjmdlTenlNiliMi,!™,) 

1 m a Mim d 

x (N*I d M' d , n d K dJd m d \T en \ NJM, pm s )* , (5) 

which is obtained from Eq. © by taking the trace over all unobserved quantum numbers of the 
electron. 

As seen from Eq. ©, the information about the states of the excited nucleus produced 
by the electron capture into the ion is now contained in the transition matrix elements 
(N2I d M d , n d K, d j d m d \T en \ NJiMi, pm s ). These matrix elements contain the wave function 
\pm s ) of a free electron with a definite asymptotic momentum. For further simplification of the 
intermediate nuclear spin-density matrix, it is therefore necessary to decompose this continuum 
wave into partial waves \eKjm), in order to apply later the standard techniques from the theory of 
angular momentum. As discussed previously lllql . however, special care has to be taken about the 
choice of the quantization axis since this directly influences the particular form of the partial wave 
decomposition. Using, for example, the direction of the electron momentum p as the quantization 



axis, the full expansion of the continuum wave function is given by UJ$\ 



\pm s ) = ^Ve iAK ^An{2l + 1) (10 l/2m s | jm s ) \eKjm s ) , (6) 



where the summation runs over all partial waves k = ±1, ±2, . . . , along all values of Dirac's angu- 
lar momentum quantum number k = ±(j+l/2) fori = j±l/2. The symbol (jim\ | j'12^12) 
generally represents the Clebsch-Gordan coefficients for the coupling of two angular momenta ji 
and j 2 to ji2- In our notation, the orbital momentum / represents the parity (—1)' of the partial 
waves \eK,jm s ), and A K is the Coulomb phase shift given by lllql 

1 f-K + iv/W\ , w ... nU + l-s) 

A K = - arg — !— - arg r (s + iv)) + '- , (7) 

2 \ s + iv I 2 



with W = Ea 2 , v = aZW/y/W 2 — 1, s = a/ k? — (aZ) 2 . Here, a is the fine- structure constant 
and E is the total electron energy. In the case of capture into ions with an initially closed shell - 
i.e., He-like - configuration, the phases can be approximated by using an effective nuclear charge 
of Z eff = Z — N b in Eq. ©, with N b being the number of bound electrons in the ion. The 
sufficiency of this approximation is confirmed by calculating the electrostatic potential induced 
by the screening electrons in the Dirac-Fock approximation II 1711 and numerically determining the 
phases for the combined nuclear and screening potentials. 

Using the decomposition © of the continuum wave function, the intermediate nuclear density 
matrix © can be rewritten in the form 

(N* d I d M d \p7 d \N* d I d M' d ) = E^' ei(A ^ AK ' ) V / (2/+l)(2/' + l) (8) 

m s Mimd kk' 

x (10 l/2m s I jm s ) (N d I d M d , n d K d j d m d \T en \ NJiMi, EKjm 3 ) 

x (I'O l/2m s I j'm s ) (N d I d M' d , n d K d j d m d \T en \ NJiMi, en'j'm s )* . 

Equation © represents the most general form of the intermediate nuclear density matrix which 
can now be used to study the properties of the excited nucleus \N d I d ) following the capture of a 
free electron. For the analysis of the radiative de-excitation of such a nucleus, however, it is more 
convenient to represent its intermediate state in terms of the statistical tensors of rank k 

Pk q (N* d I d ) = J2 (-1)^ {hM d I d - M' d I kq) (N* d I d M d \pT d \ N* d I d M f d ) . (9) 

M d M' d 

Although both the density matrix ® and the statistical tensors © contain the same physical 
information, the latter form enables one to exploit the rotational symmetry of free atoms and ions. 
By inserting the density matrix © into the definition ©, we finally obtain the statistical tensors 
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of the intermediate state as 

m s Miin d kk' M d M' d 

x (10 l/2m s | jm s ) (I'O l/2m s \ j'm s ) (I d M d I d - M' d \ kq) 



x (N d I d M d , n d K d j d m d \T en \ NJiMi, eKjm s ) 

x (N*I d M' d , n d K d3d m d \T en \ tylM, e^fm s )* . (10) 

Further evaluation of these statistical tensors within the framework of the nuclear collective 
model is discussed in the next section. The spin state of the excited nucleus is described by the 
reduced statistical tensors or alignment parameters 

which are directly related to the cross section for the population of the different nuclear magnetic 
sublevels \N*I d M d ). 

B. Calculation of the transition amplitudes 

The expression of the statistical tensors involves the matrix element of the transition operator 
T en that describes the electron-nucleus interaction in the case of electric transitions of the nucleus. 
We adopt in the following the Coulomb gauge as it allows the separation of the dominant Coulomb 
interaction between the electronic and nuclear degrees of freedom. The transition operator corre- 
sponds in lowest order to the interaction Hamiltonian H en applied in [61]: 

T en = [ d 3 r n J^4- . (12) 
J \r e -r n \ 

In the above equation, p n (r n ) is the nuclear charge density, r n denotes the nuclear coordinate and 
r e the electronic coordinate. The integration is performed over the whole nuclear volume. The 
matrix element of the transition operator T en enters the expression of the NEEC rate for electric 



transitions 



A] 



" n 



/ 



2(2Zj + 1) Mims M d m d 

(Kip \(N d I d M d , n d K d j d m d \T en \ NJiMi, pm s )\ 2 p { , (13) 



where p, is the density of the initial electronic states and J £l p denotes integration over all possible 
directions of the incoming electron. 



In order to describe the nuclear transition we use a phenomenological collective model H18H 
that interprets the characteristic band structures in the energy range up to 2 MeV in the case 
of deformed even-even nuclei as vibrations and rotations of the nuclear surface. The even-even 
nuclei have usually a low-lying 2 + first excited state, which is characterized by a strong electric 
E2 transition to the ground state. The nuclear surface is parametrized as 

oo £ 

R{6 n , <p n ) =R (l+J2 a trtXtm{Bn, ifi n )j , (14) 

1=0 m=-t 

where the amplitudes ot£ m describe the deviations of the nuclear surface with respect to the sphere 
of radius R and thus serve as collective coordinates. Using this parametrization and requiring that 
the charge be homogeneously distributed, the nuclear charge density can be written as 

Pn(r n ) = Pq® (R(9 n , <fn) - r n ) , (15) 

with the constant charge density of the undeformed nucleus given by po = if^s- Performing a 
Taylor expansion of the Heaviside function Q(R — r n ) around R we obtain 

Pn(r n ) = Po@(^o - r n ) + p 5(r n - Rq)Rq ^ a* em Y em (9 n , ip n ) + ... . (16) 

lm 

As the vibration amplitudes of the nuclear surface are supposed to be small, we neglect the terms 
of higher order in the collective coordinates a^ m . While the first term p Q(R — r n ) in the above 
equation corresponds to a round nucleus in its ground state, the second term is characterizing the 
nuclear excitation and enters the expression of the transition operator T en in Eq. (fT2l) [6]. We can 
therefore write the transition operator as 

T en = p R ya £m / drr n F — . (17) 

lm J 

It is more convenient to express the collective coordinates a £m in terms of the spherical compo- 
nents of the electric multipole transition moment Qe m , defined as tl9J] 

Qim = J tfrr*Y lm (6,<p)p n {?) . (18) 

The interaction that accounts for the electric transitions of the nucleus then yields 

T _sr^Qim f ^ 5{r n - R )Y e * m {8 n ,<p n ) 



tin 



In the calculation of the matrix element, we use the multipole expansion 



IOO L . T 

I e n| L=0 M=-L > 

where r< and r> stand for the smaller and the larger of the two radii r e and r n , respectively. 
The integration over the nuclear angular coordinates brings us to the following expression for the 
Coulomb interaction: 



« -Eo^T^( 9 «^ dry^ 5(r n - R,) . (21) 

LM J ° r > 

The matrix element of the transition operator reads 

(NjI d M d , n d K d j d m d \T en \ NJiM u 6Kjm s ) = 
lk E ^TT WW* \Qu*\ NiliMi) 

^0 LM lL + 1 

X \n d K d j d m d Yl M (9 e ,ip e ) J dr n r 2 n —^-yb{r n — Rq) 



eKjm s ) . (22) 



We write the matrix element of the electron-nucleus interaction operator as a product of the nu- 
clear and electronic parts. It is more convenient to use the reduced matrix element of the electric 



multipole operator Qlm, defined as [20] 



{N* d I d M d I Qlm \ NJ.M,) = K [ (I d M d I, - M t \ LM) 

V 2L + 1 

x (N2I d \\Q L \\NiIi) . (23) 

The modulus square of the reduced matrix element of the electric multipole operator can be related 
to the reduced electric (E) transition probability of a certain multipolarity L, 

B(EL, Ii -> I d ) = ^L-^KN^WQlWNM 2 , (24) 

whose value can be taken from experimental results. For a given multipolarity L, the matrix 
element can be written as 

{N d I d M d , n d K d j d m d \T en \ NJiMi, 6Kjm s ) = 

L 

(-lY^ +L ^ +m ^R (L+2) R L ^ K (N* d I d \\Q L \\ Nili) 

H=—L 



x ^ + 1 y (2^+1)3 fr-MihM^Lp) 

x (j - m s Jd m d \L-fj) ( Jd l/2 L0 | jl/2) , (25) 
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with the electronic radial integral Rl,k a ,k given by 

RL,K d ,n = -^tzt j£ dr e r^ +2 ^f ndKd (r e )f tK (r e ) + £n dKd (r e )# eK (r e )^ 

+R L +2 jH dr e r- e L+1 (f ndKd (r e )f eK (r e ) + 9n dKd (r e )g eK (r e )^ . (26) 

In the electronic radial integrals, g eK and f eK are the large and small radial components of the 
relativistic partial continuum electron wave function 

/ g eK {r e )tt™ s (6 e ,(p e ) \ 
(t>eKjm a {r e ) = (r e \enjm s ) = ' , (27) 

\ifen{Te)^ S K {Oe^e) ) 

and g ndKd and / ndKd are the components of the bound Dirac wave functions 

f-\ /-i ■ \ ( 9n dKd (r e )^(6 e , Ve ) \ 
4>n d K d j d m d {r e ) = (r e \n d K dJd m d } = , (28) 

\ifn d K d {r e )tt_ d Rd {e e ^ e ) J 

with the spherical spinor functions The radial integral RL,n d ,K is calculated numerically. For 
the particular case of the + — > 2 + E2 transitions, the transition amplitude reads 



(N*2M d , n d K d]d m d \T en \ 6Kjm s ) = -^==R^(-l) M <> +m ° +3jd y f 2j d ~+l 

v 125 

x (N*2 \\Q 2 \\ N.O) (j - m s Jd m d | 2 - M d ) (j d l/2 20 | jl/2) R 2)Kd}K . (29) 

The Clebsch-Gordan coefficient (j — m s j d m d \ 2 — M d ) imposes that M d = m s — m d , therefore 
M d = M' d and q = in Eq. ([T0|>. 



C. Radiative de-excitation of the nucleus 



By making use of the multipole transition amplitudes (1251) . we are able to calculate now the 
alignment parameters (fTTT) of the nucleus excited by electron capture. A detailed knowledge of 
the alignment parameters is required for the analysis of the subsequent de-excitation of the nu- 
cleus which may result in the emission of one (or several) photons until the nuclear ground state 
is reached. This photon emission is characterized (apart from its known energy) by its angular 
distribution and polarization. The relations of both of these properties to the alignment Ak of the 
excited nuclear state are well known since the early 1960s and have been discussed in detail else- 



where [21, 



22j, 



23, 



2411 . In the present work, therefore, we will restrict ourselves to a rather short 
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account of the basic formulas. For instance, the angular distribution of the gamma rays emitted in 
the transition from the exited state \N^I d M d ) to the nuclear ground state \NfIfMf) is given by 

= ^£ ( 1 + E fk(N* d I d ,N f I f )A k (N:i d )P k (co S e)) , (30) 

71 \ fc=2,4,... / 

where (Xneec is the total cross section for NEEC followed by the radiative decay of the excited 
nucleus and 9 denotes the angle of the photons with respect to the momentum p of the incoming 
electrons (chosen as the ,2-axis). As seen from Eq. d30l) . the angular dependence of the photon 
emission results from the Legendre polynomials Pk(cos9) which are weighted by the alignment 
parameters Ak(N^I d ) and the coefficients fk(N d I d , Nflf), In contrast to the alignment param- 
eters, these geometrical coefficients are independent on the nuclear excitation process and just 



account for the initial and the final nuclear states [23, 



24, 



25h, 



h(N* d I d ,N f I f ) = W±±± <^-W ( _i)W*+i (31) 



LL'pp' 



X 



V / (2L + 1)(2L' + 1) (LI L' — l\ kO) 

L V k 



■1) 



L+p+L'+p'-k 



X 



Id h If 

(N*I d ||T nr (L,p)|| N f I f y (N*I d ||r rar (L',p')|| iV^ 

Y,L P \WA\Tnr{L,p)\\N f I f )\ 2 



Here, (N d I d \\T nr (L,p)\\ Nflf) denotes the reduced matrix element for the \N^I d M d ) — > 
\NfIfMf) decay under the emission of a photon with angular momentum (or multiplicity) L and 
parity (— \) P+L . The form of the nucleus-radiation field multipole interaction operator T nr (L, p) is 
given in, e.g., Refs. flQ. 

Equations (|30l ) and (|3Tb display the general form of the photon angular distribution for the 



decay of an aligned system H23L 12511. This angular distribution includes all multipoles of the photon 
field which are allowed for the given radiative transition. In the particular case of the 2 + — > + 
nuclear decay, only the electric quadrupole transition E2 is allowed by the selection rules for 
which Eq. (|30l) simplifies to 

da^c = ^pc w{9) (32) 

dVl 4:71 

with the angular distribution given by 



W(9) ( 1 - ^-A 2 P 2 (cos6) - ^j^A 4 P 4 (cos9) ) . (33) 
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As seen from the expression above, the angular distribution of the electric quadrupole transition is 
entirely determined by the alignment parameters A\ and A4 of the excited 2 + nuclear state and the 
geometrical coefficients. In the next section we calculate these parameters and, hence, the angular 
distribution of the de-excitation photons for the electron capture into the lsi/2 and 2si/ 2 states of 
initially bare or He-like ^ 4 Yb, ™Er, *f Gd and *f Dy ions. 



III. RESULTS AND DISCUSSION 



In this section we present the alignment parameters and the angular distribution of the emitted 
photons that follow NEEC into the Is orbitals of bare ions and the 2 s states of initially He-like 
ions. We consider the even-even nuclei 7Q 4 Yb, ^Er, g| 4 Gd and <ypDy for which NEEC total cross 
sections for the capture into the K-shell have been presented in Ref. [6]. The reduced transition 
probability B(E2, 0^2) for these nuclei as well as the energies of the nuclear transitions are 
taken from Ref. [26]. The calculation of the statistical tensors involves the numerical integration 
of R,L,K d ,K- For the continuum electron we use relativistic Coulomb-Dirac wave functions, applying 
the approximation that the nucleus is a point-like charge. We assume that the free electron, which 
is far away from the ion, is not sensitive to the internal structure or size of the nucleus. The 
radial wave functions have been calculated using the same computer routine as in Ref. |]6j], and 
cross-checked with the program described in Ref. 112711 . The Coulomb phases calculated according 
to Eq. © do not include the effect of the finite nuclear size, which is expected to be negligible. 
For the continuum electrons recombining into the 2s orbital of an initially He-like ion, we use an 
effective nuclear charge number of Z cS = Z — 2. This approximation is assumed to be sufficient 



for the present leve 



of accuracy. We consider relativistic wave functions calculated with the 



GRASP92 package [28] for the bound state. The finite size of the nucleus, i.e., its non-zero radius 



R , is considered in the bound wave functions and has a sensitive effect on the inner- shell level 
energies of the bound electron. The nuclear radius in calculated according to the semi-empirical 
formula U 

R = (1.0793 A 1/3 + 0.73587) fm , (34) 

where A is the atomic mass number. 

In the case of recombination into the 2s orbital of an initially He-like ion, the phase shifts of the 
partial continuum wave functions were calculated by considering the Dirac-Fock (DF) approxima- 
tion of the Is 2 ground state seen by the free electron. As shown in Table H the difference between 
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TABLE I: Partial wave phase shifts (A K , in radian) for capture into the Is and 2s orbital of bare ions 
(Coulomb phases, third and fourth column) and into the ls 2 2s state of initially He-like ions, in different 
approximations. In the fifth column, phases for He-like ions with an effective nuclear charge Z c ^ = Z — 2 
are given and the last column (DF) contains phases corrected for bound electron screening in the Dirac-Fock 
approximation. 



Is 2s 





K 


A K (Z cS = Z) 


A K (Z cS = Z) A K (Z eff = Z-2) 


A K (DF) 


64 ^ a 


2 


2.313 


2.507 


2.524 


2.529 




-3 


-0.913 


-0.731 


-0.709 


-0.706 




2 


0.747 


2.246 


2.270 


2.273 




-3 


-2.461 


-0.983 


-0.954 


-0.952 


68° Er 


2 


0.805 


2.254 


2.277 


2.280 




-3 


-2.408 


-0.980 


-0.952 


-0.951 


^Yb 


2 


5.293 


2.205 


2.229 


2.231 




-3 


2.079 


-1.033 


-1.004 


-1.003 



such a calculation and the Z c g = Z — 2 full screening approximation is very small in the Coulomb 
phase, less than 0.01 rad and has a negligible effect on the alignment parameters. In Table U we 
present the Coulomb phases for the capture into the Is orbital and 2s orbitals of the considered 
ions, the latter calculated using both approximations. 

The values of the alignment parameters A2 and A4 of the 2 + excited nuclear states are presented 
in Table HO Capture into the Is as well as the 2s orbitals are considered. The alignment of the 
excited nuclear state characterized by these parameters gives in the second step of NEEC the 
angular distribution of the emitted radiation. In Figure |2] we present the angular distribution W{9) 
given in Eq. (1331 for the capture into the Is and 2s shells of 64 4 Gd, j^Dy, 6g°Er and 7Q 4 Yb, 
respectively. The angular patterns are similar for all the four ions, as they all involve E2 transitions 
of nuclei with near atomic and mass numbers. Both radiations following the capture into the Is 
and 2s orbitals of the ion present maxima at 9 = 45° and 9 = 135°, a pattern which significantly 
differs from that of RR [30]. While for the capture into the Is orbitals the radiation intensity 
emitted at 9 = 0°, 90°, 180° appears to be negligible, the pattern for the capture into 2s displays 
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TABLE II: Alignment parameters for the 2+ excited nuclear state formed by NEEC. We denote the capture 
orbital by nlj and E c is the energy of the continuum electron at the resonance. 





^ c (keV) 


nlj 


^2 


M 


gl^Gd 


64.005 


l s l/2 


-1.183 


1.547 


DO J 


10.318 


1/ A 


-1.185 


1.557 


170 R 
68 - Cjr 


11.350 


ls l/2 


-1.183 


1.548 


174 Yh 


4.897 


lSi/2 


-1.181 


1.542 




108.847 


2si/ 2 


-1.146 


1.383 


6 I 4 Dy 


58.164 


2si/ 2 


-1.090 


1.135 


68 7 °Er 


62.317 


2 «l/2 


-1.073 


1.056 


ro 74 Yb 


59.106 


2S1/2 


-1.029 


0.861 



larger minimum values at these angles. In contrast to NEEC, the radiative recombination of the 
free electron is dominated by the El transition and has an angular distribution (roughly given by 
a sm 2 9 function) with a maximum at 9 = 78° in the center-of-mass frame. 

As the RR and NEEC angular distributions of the emitted photons have maxima at different 
values of 6, we calculate the ratio between the two angular differential cross sections at different 
emission angles, 

^Mt'W^HL . <35) 

^ ' f — "max 

for the case of electron capture into bare ytterbium. We consider here the angles 6 max = 
45°, 78°, 135° that correspond to the maxima of NEEC and RR radiation angular distributions. 
The NEEC total cross section is convoluted with the energy distribution of the continuum elec- 
trons assuming a Gaussian width parameter of 0.1 eV. The RR angular differential cross section 
is calculated within the framework of Dirac's relativistic equation and by taking into account the 



311]. We 



higher (non-dipole) terms in the expansion of the electron-photon interaction Jlfl I30I 
assume that the RR and NEEC alignment parameters are constant on the studied energy interval 
of approximately 1 eV. 

We envisage the scenario of a possible NEEC experiment in a storage ring, in which the radi- 
ation is emitted by the nucleus of the Yb 69+ ion moving relativistically with respect to the labo- 
ratory frame. A Lorentz transformation of the NEEC and RR angular differential cross sections 
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FIG. 2: Angular distribution W{9) of photons emitted in the radiative nuclear decay of the 2 + excited state 
following NEEC for the elements g4 4 Gd, QQ 4 Dy, 68°E r > and 7Q 4 Yb, respectively. The cases of recombina- 
tion into the Is state of initially bare ions (solid lines) and into the 2s orbital of initially He-like ions (dashed 
lines) are presented. 



in the center-of-mass frame is therefore required in order to obtain the quantities in the laboratory 
system. The angular differential cross section in the laboratory system can be written as Jl^ll 

da(9) 1 da 



HOT) 

(36) 



dQ " 7 2 (1 -(3cos9) 2 dQ' ' 

where da' \6') / 'dQ! is the differential cross section in the center-of-mass system, denoted until now 
by the unprimed symbols. In our case the reduced velocity (3 is 0.138 and the Lorentz factor is 
7 = 1.009. The angle of the photons in the laboratory frame 9 is related to the one in the ion-fixed 
frame 9' by 

, cos 9 - (3 

cos 9 = . (37) 

1 — p cos 9 

As the system possesses azimuthal symmetry, cp' = (p. The angular distribution of the photons 
with respect to the laboratory system for the case of ytterbium is presented in Fig. [3] 
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FIG. 3: Angular distribution of the photons with respect to the laboratory (lab) and center-of-mass (cm) 
systems for the case of NEEC into the Is orbital of Yb 70+ . 

In Fig. |4] we present the ratio in Eq. (1351) as a function of the continuum electron energy for 
the three values of the photon emission angle 9 in the laboratory frame, for which <icx NEBC 
or da RR (9) / dVt have a maximum. In the laboratory frame, the NEEC angular distribution has 
maxima at 9 = 52° and 142°, while in the case of RR the peak is at approximately 9 = 86°. The 
ratio of the NEEC and RR angular differential cross sections is more than one order of magnitude 
larger for 9 = 52° and 142° than in the case of 9 = 86°. If the photons emitted perpendicular to the 
direction of the incoming electron are measured in an experiment, it is most likely that only the RR 
background will be detected, as the NEEC differential cross section presents a minimum at 9 = 
98°. For the emission angles 9 = 52° and 9 = 142° the RR contribution at the resonance electron 
energy E c is considerably lower in comparison to the total cross sections of the recombination 
process, 



As the ratio of the NEEC and RR angular differential cross sections is rather small, the ex- 
perimental observation of the NEEC signature is challenging. Nevertheless, knowing the angular 
pattern of NEEC is important as it provides a means of suppressing the RR background. Storage 
ring experiments focused on detecting the photons emitted in photo-recombination have the best 




(38) 



(39) 
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FIG. 4: The NEEC and RR angular differential cross sections ratio for the case of 7Q 4 Yb as a function of 
the continuum electron energy for three different photon emission angles. The NEEC total cross section 
was convoluted with a Gaussian electron energy distribution with the width parameter of 0. 1 eV for a better 
visual representation of the data. 

chances to observe the NEEC resonance at an angle of 9 = 142° for the case of E2 transitions of 
the 7Q 4 Yb nucleus. 

IV. SUMMARY 

In this work we consider the alignment of nuclear states formed in the process of nuclear 
excitation by electron capture. We calculate alignment parameters and geometric factors which 
determine the angular distribution of the subsequently emitted gamma photons. Such distribution 
functions are presented for a range of heavy elements with E2 nuclear de-excitation transitions. 
As the emission pattern of nuclear gamma photons is found to be substantially different from the 
emission characteristics of the background process of radiative recombination, our findings may 
help experimental attempts to discern NEEC from the competing process. E.g., in the case of 
the yg 4 Yb nucleus with an E2 transition, the best chance to observe the NEEC resonance is at an 
angle of 9 = 142°. This emphasizes the importance of measuring photon angle-resolved cross 
sections. Furthermore, the knowledge of the maxima and minima of the angular distribution of 
gamma radiation following NEEC is of general interest for any measurement aiming to observe 
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NEEC. 
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